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We revisit the model of the two-component plasma in a gravitational 
field, which mimics charged colloidal suspensions. We concentrate on 
the computation of the grand potential of the system. Also, a special 
sum rule for this model is presented. 



I. INTRODUCTION 

In a recent paper,0 the author presented a particular solvable model of 
Coulomb system inspired from the problem of the sedimentation equilibrium 
of charged colloidal suspensions. Although the model is much simpler than 
real colloidal suspensions, it features several properties observed in numerical 
results for more realistic mo dels! and shows some new interesting results. 

This Communication is a complement to Ref. p] in which only the density 
profiles where computed. Here we concentrate on thermodynamic quantities 
such as the grand potential. In section [TI| we briefly describe the model, in 
section |T| we compute the grand potential of the system, and in the last 
section we present a special sum rule for this model. 
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II. THE MODEL 



We consider a two-dimensional two-component plasma composed of two 
species of particles with charges ±q and masses M±. We choose a system 
of Cartesian coordinates (Oxy) in which the gravitational field is g = — gy. 
The particles are in a container of height h and infinite width. Let us define 
the inverse gravitational lengths of the particles k± = j3M±g, k = {k + + 
fc_)/2, and Sk = k + — k_, where (3 is the inverse temperature. In the grand 
canonical ensemble we define position-dependent rescaled fugacities m(r) = 
mo exp(— k±y) to account for the interaction of the particles with the external 
gravitational field. The screening length (at zero altitude) is given by 1/moP ° 

The two-component plasma is equivalent to a free Dirac field! SB when 
j3q 2 = 2. The grand partition function can be written as H = det(l + K), 
where 



K = ^+(r)-^ + m _( r )^-^j r 1 , (2-1) 

with a x , <jy, a z the Pauli matrices and (j) = a x d x + o y d y . To compute S, one 
must solve the eigenvalue problem 

KV = X9f, (2.2) 

where ^ = (ip,x) and A are the eigenvectors and eigenvalues of K. Then the 
grand potential is given by 

n = -k B rJ2H 1 + ^)- (2-3) 

A 

When an external field is acting differently over the positive and negative 
particles, as in the present case, it is useful to write the fugacities as ra±(r) = 
m(r) exp[— (±)2V(r)] where m(r) = m exp[— koy] and V(r) = Sky/4. Let us 
define the auxiliary eigenfunctions 

+ (r) = e^V(r), (2.4) 
0_(r) = e-^x(r), (2.5) 
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and the operators 



and 



a = d x + id y + d x V{r) + id y V{v) 



-d x + idy + d x V(r) - id y V(r) . 



Then the eigenvalue problem ( |2.2|) is equivalent to 



4>+(r) + A 2 a*m(r) 1 t 



mir 



,-i 



Aa 



m(r) V+( r ) 



(2.6) 
(2.7) 

(2.8) 
(2.9) 



and the boundary conditions that ip (x) on the boundary is equal to a func- 
tion which is analytic (anti-analytic) outside the container and vanishes at 
infinity.BB This method is general and can be applied to others models of 
two-component plasma in an external field. 



III. THE GRAND POTENTIAL 



In the present case, we look for solutions of (2.8) of the form 0+(r) 
h(y) exp(ikx — k^y). From (|2.8| ) we find the equation for h(y) 



(k + 5k/A)k - (jfc + 5k/ A) 2 - \- 2 m 2 e- 2koy ] h(y) + k h'(y) + h"(y) = , 

(3.1) 

which is very similar to the equation (3.2) of Ref. |l| satisfied by the Green 
function g ++ . The solution is 



%) 



AI U 



-koy s 



\k c 



+ BK U 



m e 



\k 



-hoy/2 



(3.2) 



with v — \k — (k-/2)\/ko, and A and B constants of integration. 

In the present case, the boundary conditions become <ft + (x,y = 0) = 
and 0_ (x, y — h) — if k > 0, and (p + (x,y = h) = and <fi-(x,y = 0) = if 
k < 0. These conditions give a system of two linear homogeneous equations 
for the constants A and B. Writing that the discriminant of this system 
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must be zero, in order to have non-trivial solutions, gives the equation that 
determines A. Let a be the sign of k — (&_/2). If k > 0, A is a solution of 



m 



moe 



-koh ' 



" \Xk n J \ Xk Q 



and if k < 0, A is given by 



moe 



-koh s 



Xkn 



m e 



-koh " 



\k n 



K v - X 



/mo 
\Xk n 



1^ K, 1 m ° 6 



\Xkr 



m 
Xk 



-koh s 







(3.3) 



Xk n 







(3.4) 



These equations have an infinite number of solutions. For each k, we may 
index the different values of A by an integer £: X = A^. Then, the grand 
potential by unit length in the x-direction is 



k B T 
2vr 



J 5>(1 + A M ) dk. 



(3.5) 



The sum over £ can be done explicitly,^ i for example in the case where 
k < 0, by noting that the zeros of the entire function 



kn 



I u (m e fc °V^o)^-i("W&o) 



+I u - 1 (m z/k )K u ( y m e k ° h z/k n ) 

are 1/Xki and that /(0) = 1, so f(z) = Y[(l — Xktz). Then, 

e 

$>(1 + A M ) = hint! + A M ) = ln/(-l) . 



(3.6) 



(3.7) 



The same can be done in the case k < 0. Then, the grand potential can 
be expressed as three integrals for k in the domains ]— oo,0], [0, k_/2] and 
[k-/2, +oo[. After a change of variable in the integrals, the grand potential 
can be written in a more compact way as 



puj = ( 

2n J-k-/2k 



+oo 



In 



k+/2k 



e- koh{v+1) ^{l„(^)K v+1 (^e- k ° h ) (3.8) 



dv . 



From now on we consider the simpler case when the container has infinite 
width (h — > oo). Equation (|3.8|) becomes 



p u = -±L(l +1 )ln 

2% J-k-/2k J~k + /2k 



( 2 -^\ r(„ + i)i„(=l 

V m I \k 



dv. (3.9) 



The two-component plasma at f3q 2 = 2 has logarithmic divergences due 
to the collapse of pairs of particles of opposite sign. This can be avoided by 
introducing hard core particles of radius R. The integrals in ( |3.8| ) and ( |3.9| ) 
are divergent and must be cutoff, the upper limit of integration becomes 

^max = (koR) ■ 

We can have an asymptotic expression for the grand potential in the 
usual physical case where k <C m , using the Debye expansioni of the Bessel 
functions. The calculations are very similar to those of Ref. || for a different 
problem. One finds 

" = -inrPo + o( 1 ¥) , (3.10) 



2k \ k 

where po is the pressure of a two-component plasma without an external field. 
It is interesting to notice that this is the same grand potential of a system 
without gravity, but confined in a container of height (2ko)~ 1 . Furthermore, 
from Ref. we know that most of the particles are in fact in the region 
< y < (2A; ) _1 since for intermediate altitudes (the neutral zone) the 
density profiles decay as exp(— 2k y). 



IV. A SPECIAL SUM RULE 

We can write a sum rule for this system by computing in two different 
ways the force exerted on the bottom of the container. This force is the 
weight of the column of fluid over the base of the container. On the other 
hand, it is also given by the pressure at zero altitude which can be computed 
by means of the contact theorem: it is the density at zero altitude times fc_BT. 
This gives the sum rule 

/ r+oo r+oo \ 

P+ (0)+p-(p) = Pg\M + J o p + (y)dy + M.J Q P -(y)dy), (4.1) 
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where p± are the individual densities of the positive and negative particles 
respectly. This sum rule can be written in a simpler way; because the system 
is neutral, both integrals in ( j4.1|) are equal to n/2, where n is the total number 
of particles. Then, the sum rule becomes 

p+(0) + p_(0) = k n. (4.2) 

We expect this to be true at any temperature. We can test the sum rule 
in the present case when f3q 2 = 2. The total number of particles can be 
computed from the grand potential: 

n = -Prno-j^ ■ (4.3) 

From equation (|3.9|) one finds 



ran f +oc ^+i(f a ) 

n ='-!^( + J° dv. (4.4) 

2tT i-fc_/2fc J-k+/2k J I v (*g) V ' 

On the other hand, from equation (3.10) of Ref. [I] we can compute the den- 
sities at zero altitude: 



The sum rule (14. 21) is verified. 
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